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The Feynman, Metropolis and Teller treatment of compressed atoms is extended to the relativistic 
regimes. Each atomic configuration is confined by a Wigner-Seitz cell and is characterized by a 
positive electron Fermi energy. The non-relativistic treatment assumes a point-like nucleus and 
infinite values of the electron Fermi energy can be attained. In the relativistic treatment there 
exists a limiting configuration, reached when the Wigner-Seitz cell radius equals the radius of the 
nucleus, with a maximum value of the electron Fermi energy (Eg) max , here expressed analytically in 
O ' the ultra-relativistic approximation. The corrections given by the relativistic Thomas-Fermi-Dirac 

exchange term are also evaluated and shown to be generally small and negligible in the relativistic 
high density regime. The dependence of the relativistic electron Fermi energies by compression for 
selected nuclei are compared and contrasted to the non-relativistic ones and to the ones obtained in 
the uniform approximation. The relativistic Feynman, Metropolis, Teller approach here presented 
overcomes some difficulties in the Salpeter approximation generally adopted for compressed matter 
in physics and astrophysics. The treatment is then extrapolated to compressed nuclear matter 
cores of stellar dimensions with A ~ (rapi anc k/m n ) 3 ~ 10 57 or M core ~ Mq. A new family of 
equilibrium configurations exists for selected values of the electron Fermi energy varying in the 
■ range < Eg < (£f ) mal . Such configurations fulfill global but not local charge neutrality. They 

have electric fields on the core surface, increasing for decreasing values of the electron Fermi energy 
q , reaching values much larger than the critical value E c = m^c 3 /(eh), for Eg = 0. We compare and 

5_l . contrast our results with the ones of Thomas-Fermi model in strange stars. 

CZ2 ! 

PACS numbers: 31.15.bt, 71.10.Ca, 26.60.Dd 

Keywords: Thomas-Fermi model; Degenerate Fermi gas; Neutron star core. 

> 

CN . I. INTRODUCTION 

In a classic article Baym, Bethe and Pethick [1| presented the problem of matching, in a neutron star, a liquid 
core, composed of N n neutrons, N p protons and N e electrons, to the crust taking into account the electrodynamical 
and surface tension effects. After discussing the different aspects of the problem they concluded: The details of this 
picture requires further elaboration; this is a situation for which the Thomas- Fermi method is useful. This statement, 
• • ■ in first instance, may appear surprising: the Thomas-Fermi model has been extensively applied in atomic physics 
(see e.g Gombas @, March Q, Lundqvist and March 0), also has been applied extensively in atomic physics in its 
relativistic form (see e.g Ferreirinho, RufHni and Stella Q, Rumni and Stella as well as in the study of atoms with 
heavy nuclei in the classic works of Migdal, Popov and Voskresenskii 0, Q . Similarly there have been considerations 
of relativistic Thomas-Fermi model for quark stars pointing out the existence of critical electric fields on their surfaces 
(see e.g. Alcock, Farhi, Olinto Usov Similar results have also been obtained by Alford et al. [11] in the 

transition at very high densities, from the normal nuclear matter phase in the core to the color-flavor-locked phase of 
quark matter in the inner core of hybrid stars. No example exists to the application of the electromagnetic Thomas- 
Fermi model for neutron stars. This problem can indeed be approached with merit by studying the simplified but 
rigorous concept of a nuclear matter core of stellar dimensions which fulfills the relativistic Thomas-Fermi equation 
as discussed in [12, EH . As we will see this work leads to the prediction of the existence of a critical electric field at 
the interface between the core and the crust of a neutron star. 

In pjl, [l3| we have first generalized the treatment of heavy nuclei by enforcing self-consistently the condition of beta 
equilibrium in the relativistic Thomas-Fermi equation. Using then the existence of scaling laws we have extended the 
results from heavy nuclei to the case of nuclear matter cores of stellar dimensions. In both these treatments we had 
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there assumed the Fermi energy of the electrons E^ — 0. The aim of this article is to proceed with this dual approach 
and to consider first the case of compressed atoms and then, using the existence of scaling laws, the compressed 
nuclear matter cores of stellar dimensions with a positive value of their electron Fermi energies. 

It is well known that Salpeter has been among the first to study the behavior of matter under extremely high 
pressures by considering a Wigner-Seitz cell of radius Rws 0- Salpeter assumed as a starting point the nucleus 
point-like and a uniform distribution of electrons within a Wigner-Seitz cell. He then considered corrections to the 
above model due to the inhomogeneity of electron distribution. The first correction corresponds to the inclusion of 
the lattice energy Eq — — (9NpO)/ (lORws), which results from the point-like nucleus-electron interaction and, from 
the electron-electron interaction inside the cell of radius Rws- The second correction is given by a series-expansion 
of the electron Fermi energy about the average electron density n e given by the uniform approximation. The electron 
density is then assumed equals to n e [l + e(r)] with e(r) considered as infinitesimal. The Coulomb potential energy 
is assumed to be the one of the point-like nucleus with the uniform distribution of electrons of density n e thus 
the correction given by e(r) is neglected on the Coulomb potential. The electron distribution is then calculated at 
first-order by expanding the relativistic electron kinetic energy about its value given by the uniform approximation 
considering as infinitesimal the ratio eV/Ef between the Coulomb potential energy eV and the electron Fermi energy 
E^ — \/[cPF (r)} 2 + m 2 c' i — m e c 2 — eV. The inclusion of each additional Coulomb correction results in a decreasing 
of the pressure of the cell Ps by comparison to the uniform one (see [l5[ for details) . 

It is quite difficult to assess the self-consistency of all the recalled different approximations adopted by Salpeter. 
In order to validate and also to see the possible limits of the Salpeter approach, wc consider the relativistic gener- 
alization of the Feynman, Metropolis, Teller treatment [l6| which takes automatically and globally into account all 
electromagnetic and special relativistic contributions. We show explicitly how this new treatment leads in the case of 
atoms to electron distributions markedly different from the ones often adopted in the literature of constant electron 
density distributions. At the same time it allows to overcome some of the difficulties in current treatments. 

Similarly the point-like description of the nucleus often adopted in literature is confirmed to be unacceptable in the 
framework of a relativistic treatment. 

In Sec. |TT] we first recall the non- relativistic treatment of the compressed atom by Feynman, Metropolis and Teller. In 
Sec. IIIII we generalize that treatment to the relativistic regime by integrating the relativistic Thomas- Fermi equation, 
imposing also the condition of beta equilibrium. In Sec. IIVI we first compare the new treatment with the one 
corresponding to a uniform electron distribution often used in the literature and to the Salpeter treatment. We also 
compare and contrast the results of the relativistic and the non-relativistic treatment. 

We then proceed to analyze the case of compressed nuclear matter cores of stellar dimensions. 

In Sec. [V] using the same scaling laws adopted in 0, 0] we turn to the case of nuclear matter cores of stellar 
dimensions with mass numbers A w (?npianck/w n ) 3 ~ 10 57 or M core ~ Mq where m n is the neutron mass and 
"IPlanck = (fic/G) 1 / 2 is the Planck mass. Such a configuration present global but not local charge neutrality. Analytic 
solutions for the ultra-relativistic limit are obtained. In particular we find: 

1) explicit analytic expressions for the electrostatic field and the Coulomb potential energy, 

2) an entire range of possible Fermi energy for the electrons between zero and a maximum value (E^) max , reached 
when Rws = Rc> which can be expressed analytically, 

3) the explicit analytic expression of the ratio between the proton number N p and the mass number A when 
Rws — Rc- 

We turn then in Sec. IVII to the study of the compressional energy of the nuclear matter cores of stellar dimensions 
for selected values of the electron Fermi energy. We show that the solution with Ef = presents the largest value of 
the electrodynamical structure. 

We finally summarize the conclusions in Sec. IVIII 

II. THE THOMAS-FERMI MODEL FOR COMPRESSED ATOMS: THE 
FEYNMAN-METROPOLIS-TELLER TREATMENT 

A. The classical Thomas-Fermi model 

The Thomas-Fermi model assumes that the electrons of an atom constitute a fully degenerate gas of fermions 
confined in a spherical region by the Coulomb potential of a point-like nucleus of charge +eN p [l7l Il8j|. Feynman, 
Metropolis and Teller have shown that this model can be used to derive the equation of state of matter at high 
pressures by considering a Thomas- Fermi model confined in a Wigner-Seitz cell of radius Rws 0]- 
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We recall that the condition of equilibrium of the electrons in an atom, in the non-relativistic limit, is expressed by 

where m e is the electron mass, V is the electrostatic potential and is their Fermi energy. 
The electrostatic potential fulfills, for r > 0, the Poisson equation 

V 2 V = 47ren e , (2) 

where the electron number density n e is related to the Fermi momentum by 



3ir 2 h 3 



(3) 



For neutral atoms and ions n e vanishes at the boundary so the electron Fermi energy is, respectively, zero or negative. 
In the case of compressed atoms n e does not vanish at the boundary while the Coulomb potential energy eV is zero. 
Consequently is positive. 
Assuming 

eV(r)+EF = e 2 iV p ^, (4) 
we obtain the following expression for the electron number density 

where the new dimensionless radial coordinate r\ is given by r — br\, where 

&=(37r) V3_^ _L i (6) 



Eq. ([2]) can be then written in the form 

d 2 (/)(r)) <X?7) 3/2 



dr] 2 rj 1 / 2 



(7) 



which is the classic Thomas- Fermi equation [jq |. A first boundary condition for this equation follows from the 
point-like structure of the nucleus 

m = i. (8) 

A second boundary condition comes from the conservation of the number of electrons N e — J^ ws 47m e (r)r 2 dr 

1 - j^- = 4>(Vo) ~ Votfivo), (9) 

where t/q = Rws/b defines the radius Rws of the Wigner-Seitz cell. In the case of compressed atoms N e = N p so the 
Coulomb potential energy eV vanishes at the boundary Rws- As a result, using Eqs. ^ and ([3]), the Fermi energy 
of electrons is given by 

F _ N p e 2 4>{ m ) 

Ee ~ ~ m~' (10) 

Therefore in the classic treatment 770 can approach zero and consequently the range of the possible values of the Fermi 
energy extends from zero to infinity. 

The results are summarized in Figs. Q]and[2] 
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FIG. 1: Physically relevant solutions of the Thomas- Fermi Equation with the boundary conditions ((8| and (J5J. The curve 
1 refers to a neutral compressed atom. The curve 2 refers to a neutral free atom. The curve 3 refers to a positive ion. The 
dotted straight line is the tangent to the curve 1 at the point (rjo, 4>(t)o)) corresponding to overall charge neutrality (see Eq. 
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FIG. 2: The electron Fermi energies for iron, in units of the electron mass, are plotted as a function of the dimensionless 
compression parameter rjo- Points refer to the numerical integrations of the Thomas- Fermi equation ([7Jl performed originally 
by Feynman, Metropolis and Teller in fl6l ]. 

B. The Thomas- Fermi-Dirac model 

Dirac has introduced modifications to the original Thomas-Fermi theory to include effects of exchange [l9j . In this 
case the condition of equilibrium of the electrons in the atom is generalized as follows 



f ^2 



Pi) 

2m r - 



eV--cPf = Ei 



(11) 



where as usual a = e 2 /he denotes the fine structure constant. 

The electron number density is now connected to the Coulomb potential energy by 
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—m e c + a/ y—m e c 2 J + 2m e c 2 (eV + E^) 



(12) 



Assuming 



z \7T / r 
and r = bq, the Poisson equation can be written as 



(13) 



d 2 m 

drj 2 



V v 



1/2 



(14) 



where b is given by Eq.© and d = (3/(32tt 2 )) 1 / 3 (1/N p ) 2 / 3 . The boundary condition for Eq. CE! are 0(0) = 1 and 
W07o) = 0(?7o)- 



III. THE RELATIVISTIC GENERALIZATION OF THE FEYNMAN-METROPOLIS-TELLER 

TREATMENT 



A. The relativistic Thomas-Fermi model for atoms 



In the relativistic generalization of the Thomas-Fermi equation the point-like approximation of the nucleus must 
be abandoned [5, 6] since the relativistic equilibrium condition 

El = ^J(Pfc) 2 +m 2 c 4 - m e c 2 - eV(r) , (15) 

which generalizes the Eq. ([1}, would lead to a non-integrable expression for the electron density near the origin. 
Consequently we adopt an extended nucleus. Traditionally the radius of an extended nucleus is given by the phe- 
nomenological relation R c = tqA 1 ^ 3 where A is the number of nucleons and tq = 1.2 x 10~ 13 cm. Further it is possible 
to show from the extremization of the semi-empirical Weizsacker mass-formula that the relation between A and N p 
is given by 



N p = 



2 
A 



1 



200 A 1 / 3 



which in the limit of small A gives 



p 2 



(16) 



(17) 



In [TH we have relaxed, for E^ — 0, the condition N p w A/ 2 (adopted, for example, in Migdal, Popov and Voskresenski 
[H) as well as the condition N p = [2/A + 3/(200A 1/3 )]" 1 (adopted for example in Ferreirinho, Ruffini and Stella 
imposing explicitly the beta decay equilibrium between neutron, protons and electrons. 

In particular, following the previous treatments (see e.g. |l3j ), we have assumed a constant distribution of protons 
confined in a radius R c defined by 



R r = A- 



JVV3 
1 y p 1 



(18) 



where is the pion mass and A is a parameter such that A « 1 (A < 1) corresponds to nuclear (supranuclear) 
densities when applied to ordinary nuclei. Consequently, the proton density can be written as 



n p (r) 



InR 3 



0(R C ~ r) 



mtc 



1 



4tt h 3 A 3 



e(Rc-r), 



(19) 
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where 6{x) is the Heaviside function which by definition is given by 

6(x) = 

The electron density is given by 



0, x < 0, 

1, x > 0. 



F\3 



n e (r) = 



1 



(Pf) 

37T 2 n 3 37T 2 ft 3 C 3 



[e 2 V 2 {r) + 2m e c 2 eV{r)] 



3/2 



where V is the Coulomb potential. 

The overall Coulomb potential satisfies the Poisson equation 

V 2 F(r) = -47re [n p (r) - n e (r)} , 



(20) 



(21) 



(22) 



with the boundary conditions V(oo) = (due to global charge neutrality) and finiteness of V(0). 

By introducing the dimensionless quantities x = r/X„, x c = i? c /A w and \/ r = eV(r)/(ch) with A,r = h/im^c), and 
replacing the particle densities (p~9|) and (|26|) into the Poisson equation ([22]) we obtain the relativistic Thomas-Fermi 
equation 



1 d 2 X (x) 
3x dx 2 



4q 
9tt 



x 2 (^) 



2— - 

a; 



.3/2 



(23) 



where x(0) = 0, x(xoo) = 0. The neutron density n n (r), related to the neutron Fermi momentum = (37r 2 ?i 3 n ri ) 1 / 3 , 
is determined, as in the previous case [l3| , by imposing the condition of beta equilibrium 



K = \ (PFc) 2 +m 2 n c*-m n c 2 



= J (P/c) 2 + m^c 4 - m p c 2 + eV(r), 



(24) 



which in turn is related to the proton density n p and the electron density by Eqs. (|21j) . (|22j) . Integrating numerically 
these equations we have obtained a new generalized relation between A and N p for any value of A. In the limit of 
small A this result agrees with the phenomenological relations given by Eqs. ([TBI [T7| . as is clearly shown in Fig. ^ 
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Mass number A 



FIG. 3: The A-N p relation at nuclear density (solid line) obtained from first principles compared with the phenomenological 
expressions given by N p ~ A/2 (dashed line) and Eq. (|16|) (dotted line). The asymptotic value, for A (mpi anc k/m„) 3 , is 
N p » 0.0046A. 
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B. The relativistic Thomas- Fermi model for compressed atoms 



We turn now to the case of compressed atoms in which the electron Fermi energy is positive. The relativistic 
generalization of the equilibrium condition ([1} now reads 



^(Pfc) 2 +m 2 e c 4 - m e c 2 - eV(r) > . 



(25) 



Adopting an extended-nucleus with a radius given by Eq. (|18j) and a proton density given by Eq. (|19j) the Poisson 
equation (j2"2"|) , with the following electron density 



(Pf) 3 

3ir 2 h 3 3vr 2 /i 3 c 3 



n e {r) 



1 



e 2 V 2 {r) + 2m e c 2 eV(r) 



3/2 



(26) 



can be written as 



1 d 2 X (x) 
3x dx 2 



A 3 



9(x c 



4a 
9^ 



-, 3/2 



(27) 



where x = r/X n , x c — Rc/X^, \jr = eV(r)/(cti), = K/{m„c) and eV — eV + E[ . The equation (|2"T|) has to be 
integrated with the boundary conditions x(0) = 0, x( x ws) = %wsx'( x ws)) x\vs = i?irsA,. 

The neutron density n„(r), related to the neutron Fermi momentum P^ — (37r 2 ft 3 n„) 1 / 3 , is determined by imposing 
the condition of beta equilibrium 



El = iJ(P r ?c) 2 + m 2 n c±-m n c 2 

= tJ(PFc) 2 + m 2 c 4 - m p c 2 + eV(r) + E\ 



(28) 



Using this approach, it is then possible to determine the beta equilibrium nuclide as a function of the density of 
the system. Infact, as suggested by Hund [2(| and Landau [2l|, when the electron Fermi energy is sufficiently high, 
electrons can be absorbed by protons and converted to neutrons in inverse beta decay p + e~ — > n + v e because 

iripC 2 + eV(r) + E^ holds. The condition of equilibrium (|28]l is crucial, for 



the condition E% < ^(Pfc) 2 + m 2 c 4 

example, in the construction of a self-consistent equation of state of high energy density matter present in white 
dwarfs and neutron star crusts |15| . In the case of zero electron Fermi energy the generalized A — N p relation of Fig. 
([3]) is obtained. 



C. The relativistic Thomas-Fermi-Dirac model for compressed atoms 



We now take into account the exchange corrections to the relativistic Thomas- Fermi equation (|27j) . In this case we 
have (see [§[ for instance) 



El = J (cPf) 2 + mlc 4 - m e c 2 - eV cP/ = constant . 



Introducing the function x( r ) as before 



El +eV = eV= Hc- 



(29) 



(30) 



we obtain the electron number density 



e 37T 2 ft 3 C 

where 7 = {a/ii)/{l — a 2 /n 2 ) 



^3 <| 7 (m e c 2 + eV^j + (eV^j + 2m e c 2 eV 



1/2 



(1 + 7 2 )(m e c 2 + eV) 2 - m^c 4 
(m e c 2 + eV) 2 — m 2 c 4 



1/2' 



(31) 
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If we take the approximation 1 + -f 2 1 the above equation becomes 

1 



3ir 2 h s 



- I 7 (m e c 2 + eV^j + (eV) + 



2m e c eV 



1/2 ^ 



(32) 



The second term on the right-hand-side of Eq. (|32[) has the same form of the electron density given by the rela- 
tivistic Thomas-Fermi approach without the exchange correction (|26p and therefore the first term shows the explicit 
contribution of the exchange term to the electron density. 

Using the full expression of the electron density given by Eq. (|3Tj) we obtain the relativistic Thomas-Fermi-Dirac 
equation 



1 d 2 x(x) a M _ 

3x dx 2 ~ A3 U[Xc X) 




( m e 








m e v 
12—- 


1/2 ■ 










TOtt X 





(1 + 7 2 )(m e /m ;r + x/x) 2 - (rrie/m^) 2 



(me/rriTr + xl x ) 2 ~ {me/m^) 2 



1/2- 



(33) 



which by applying the approximation 1 
1 d 2 X (x) 



1 becomes 



3a; dx 2 



4a J (m e y n 

97T V Win X , 





[(f) 2 - 


2 m e X 


1/2 ] 






m-n X 





(34) 



The boundary conditions for Eq. (f3"3"f are %(0) = and x( x ws) — x wsx'( x ws)- The neutron density can be 
obtained as before by using the beta equilibrium condition (|28p with the electron Fermi energy given by Eq. (I29[) . 

In Fig. H] we show the results of the numerical integration of the relativistic Thomas- Fermi equation (|2"T|) and of 
the relativistic Thomas-Fermi-Dirac equation (|33[) for helium, carbon and iron. In particular, we show the electron 

Fermi energy multiplied by N p 4 ^ 3 as a function of the ratio Rws/Rc between the Wigner-Seitz cell radius Rws arid 
the nucleus radius R c given by Eq. (p~8|) . 

The effects of the exchange term are appreciable only in the low density (low compression) region, i.e. when 
Rws » R c (see Fig. [4j). We can then conclude in total generality that the correction given by the Thomas-Fermi- 
Dirac exchange term is, small in the non-relativistic low compression (low density) regime, and negligible in the 
relativistic high compression (high density) regime. 



IV. COMPARISON AND CONTRAST WITH APPROXIMATE TREATMENTS 



There exists in the literature a variety of semi-qualitative approximations adopted in order to describe the electron 
component of a compressed atom (see e.g. Biirvcnich, Mishustin and Greiner [22j for applications of the uniform 
approximation and e.g Chabrier and Potekhin (23|, Potekhin, Chabrier and Rogers [24|), Haensel and Zdunik [25| . 
[261 ] . [28| . Douchin and Haensel [27j for applications of the Salpeter approximate treatment). 

We shall see how the relativistic treatment of the Thomas-Fermi equation affects the current analysis of compressed 
atoms in the literature by introducing qualitative and quantitative differences which deserve attention. 



A. Relativistic FMT treatment vs. relativistic uniform approximation 



One of the most used approximations in the treatment of the electron distribution in compressed atoms is the one 
in which, for a given nuclear charge +eN p , the Wigner-Seitz cell radius Rws is defined by 

Air 

N p = —R ws n e , (35) 

where n e = (P e F ) 3 /(37r 2 S 3 ). The Eq. ([33]) ensures the global neutrality of the Wigner-Seitz cell of radius Rws 
assuming a uniform distribution of electrons inside the cell. 

We shall first compare the Feynman-Metropolis- Teller treatment, previously introduced, with the uniform approx- 
imation for the electron distribution. In view of the results of the preceding section, hereafter we shall consider the 
non-relativistic and the relativistic formulation of the Feynman-Metropolis- Teller treatment with no Thomas-Fermi- 
Dirac exchange correction. 
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FIG. 4: The electron Fermi energies in units of m^<?Np is plotted for helium, for carbon and for iron, as a function of 
the ratio Rws /Rc in the relativistic Feynman-Metropolis- Teller (FMT) treatment with and without the Thomas- Fcrmi-Dirac 
exchange effects. Here Rws denotes the Wigner-Seitz cell radius and R c is the nucleus radius as given by Eq. (|18[) . It is clear 
that the exchange terms are appreciable only in the low density region and are negligible as Rws — > Rc 



In Fig. [5] we have plotted the electron number density obtained from Eq. (|26p where the Coulomb potential is 
related to the function \i which is obtained from numerical integration of the relativistic Thomas-Fermi equation 
(|27p for different compressions for helium and iron. We have normalized the electron density to the average electron 
number density n = SN e /(AirR ws ) = 3N p /(4:TtR ws ) as given by Eq. (|35|) . 

We can see in Fig. [5] how our treatment, based on the numerical integration of the relativistic Thomas- Fermi 
equation (1271) and imposing the condition of beta equilibrium (|28[) . leads to electron density distributions markedly 
different from the constant electron density approximation. 

From Eqs. (|25|) . (|35|) and taking into account the global neutrality condition of the Wigner-Seitz cell eV(Rws) = 0, 
the electron Fermi energy in the uniform approximation can be written as 



E 




(36) 



We show in Fig. [6] the electron Fermi energy as a function of the average electron density uq = 3N e / (AnR ws ) = 
2>N p /(A-kR ws ) in units of the Bohr density neohr = 3/(4vr.Rg h ) where i?Bohr = ^ 2 /(e 2 m e ) is the Bohr radius. For 
selected compositions we show the results for the relativistic Feynman-Metropolis- Teller treatment, based on the 
numerical integration of the relativistic Thomas- Fermi equation (|27p. and for the relativistic uniform approximation. 

As clearly shown in Fig. [5] and summarized in Fig. [5] the relativistic treatment leads to results strongly dependent 
at low compression from the nuclear composition. The corresponding value of the electron Fermi energy derived from 
a uniform approximation overevaluates the true electron Fermi energy (see Fig. [5]). In the limit of high compression 
the relativistic curves asymptotically approach the uniform one (see also Fig. [SJ. 

The uniform approximation becomes exact in the limit when the electron Fermi energy acquires its maximum value 
as given by 




3tt 2 



2/?, 



Np 
A 



2/3 



(37) 



which is attained when Rws coincides with the nuclear radius R c . Here, the maximum electron Fermi energy (|37[) 
is obtained by replacing in Eq. (|3"6")l the value of the normalized Wigner-Seitz cell radius xws — x c — R c /Xtt ~ 
[(3/2)7r] 1 / 3 A 1 / 3 , where we have approximated the nuclear density as n nuc w (1/2)A~ 3 . 
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\ Uniform approximation 
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FIG. 5: The electron number density n e in units of the average electron number density no = 3A r e /(47r_R^ iS ) is plotted as a 
function of the dimensionless radial coordinate x = r/X n for the selected compressions xws = 9.7 (upper panels), xws = 3 x 10 3 
(middle panels) and xws = 10 4 (bottom panels), in both the relativistic Feynman, Metropolis, Teller approach and the uniform 
approximation respectively for Helium (panels on the left) and Iron (panels on the right). 



B. Relativistic FMT treatment vs. Salpeter approximate treatment 

Corrections to the uniform distribution were also studied by Salpeter [3] and his approximations are largely applied 
in physics (see e.g. Chabrier and Potekhin (23[, Potekhin, Chabrier and Rogers [24|]) and astrophysics (see e.g. Haensel 
and Zdunik [H[, [26|, @, Douchin and Haensel [27| ). 

Keeping the point-like nucleus assumption, Salpeter 14] studied the corrections to the above models due to the 
inhomogeneity of the electron distribution inside the Wigner-Seitz cell. He expressed an analytic formula for the 
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4 Rel. uniform approximation 

10 - Rel. uniform approximation with point-like nucleus 
Rel. FMT treatment for He 



Rel. FMT treatment for Fe 




FIG. 6: The electron Fermi energies E e in the relativistic Feynman-Metropolis- Teller (FMT) treatment and in the uniform 
approximation in units of the pion rest mass, are plotted as a function of the average electron density no = SN e /(^TvRws) m 
units of the Bohr density riBohr = 3/(47r.R Bohr ) where i?Bohr = h 2 /(e 2 m e ) is the Bohr radius. The filled circles correspond to 
the case of a relativistic uniform approximation with a point-like nucleus. In such a case the electron Fermi energy can reach 
arbitrary large values as Rws — ► 0. The arrow indicates the value of the maximum electron Fermi energy as given by Eq. (I37[l . 



total energy of a Wigner-Seitz cell based on Coulomb corrections to the uniform distribution of electrons. The first 
correction corresponds to the inclusion of the lattice energy Eq = —(SN^a)/ (lORws), which results from the point- 
like nucleus-electron interaction and, from the electron-electron interaction inside the cell of radius R ws - The second 
correction is given by a series-expansion of the electron Fermi energy about the average electron density n e given by 
Eq. (j3"51 the uniform approximation n e = 3Z/(47ri?^ s ). The electron density is then assumed equals to n e [l + e(r)] 
with e(r) considered as infinitesimal. The Coulomb potential energy is assumed to be the one of the point-like nucleus 
with the uniform distribution of electrons of density n e given by, thus the correction given by e(r) is neglected on the 
Coulomb potential. The electron distribution is then calculated at first-order by expanding the relativistic electron 
kinetic energy 



e fe = yj [cP/»] 2 + m 2 e c A - m e c 2 



= \J (37r 2 n e ) 2 /3[i + e ( r )]2/3 + m 2 c 4 „ TOeC 2 ; 
about its value given by the uniform approximation 



(37r 2 n< 



N2/3 



m e c 



(38) 



(39) 



considering as infinitesimal the ratio eV/E^ between the Coulomb potential energy eV and the electron Fermi energy 
JSf = y/[cP[ (r)] 2 + m 2 c 4 - m e c 2 - eV. 

The effect of the Dirac electron-exchange correction on the equation of state was also considered by Salpeter 
[ljj . However, adopting the general approach of Migdal et al. Q, these effects are negligible in the relativistic regime 
(see Subsec.Hirni). 

The inclusion of each additional Coulomb correction results in a decreasing of the pressure of the cell Ps (see [HI 
for details). However, despite to be very interesting in identifying piecewise contributions to the total pressure, the 
validity of the Salpeter approach needs a verification by a more general treatment. For instance, the failure of the 
Salpeter formulas can be seen at densities of the order of ~ 10 2 — 10 3 g cm~ 3 for nuclei with large N p , as in the 
case of iron, where the pressure becomes negative (see Table (jT|). Therefore, the problem of solving the relativistic 
Thomas-Fermi equation within the Feynman, Metropolis, Teller approach becomes a necessity, since this approach 
gives all the possible Coulomb and relativistic contributions automatically and correctly. 
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TABLE I: Pressure for Iron as a function of the density p in the uniform approximation (P), in the Salpeter approximation 
(Ps) and in the relativistic Feynman-Metropolis- Teller approach (Pf MTrei ) • Here xs = P^s / ( m e c ) , XFMTrel = Pf /(m e c) are 
respectively the normalized Salpeter Fermi momentum and the relativistic Feynmann-Metropolis- Teller Fermi momentum. 
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C. Relativistic FMT treatment vs. non-relativistic FMT treatment 

In order to compare and contrast the Fermi energy of a compressed atom in the non-relativistic and the relativistic 
limit we first express the non-relativistic equations in terms of the dimensionless variables used for the relativistic 
treatment. We then have 

Ajt r cn 

and the non-relativistic limit of Eq. ([27]) becomes 

f ^ (41) 



dx 2 Sir \m n J x 1 / 2 

with the boundary conditions 

x (0) = aN p , xwsx(xws)' = x( x ws), (42) 

and dimensionless variable xws = Rws/^-k- 

In these new variables the electron Fermi energy is given by 

E! = (43) 
xws 

The two treatment, the relativistic and the non-relativistic one, can be now directly compared and contrasted by 
using the same units (see Fig. [7]). 
There are major differences: 
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FIG. 7: The electron Fermi energies in units of m 7r c 2 iVp / ' 3 for helium, for carbon and for iron are plotted as a function of the 
ratio Rws/Rc respectively in the non-relativistic and in the relativistic Feynman-Metropolis- Teller (FMT) treatment without 
the Thomas-Fermi-Dirac exchange effects. Here Rws is the radius of the Wigner-Seitz cell and R c is the radius of the nucleus 
given by Eq. (|18|) . The relativistic treatment leads to results of the electron Fermi energy strongly dependent on the nuclear 
composition and systematically smaller than the non-relativistic ones, which can attain arbitrary large values as the point-like 
nucleus is approached. 



1) The electron Fermi energy in the relativistic treatment is strongly dependent on the nuclear composition, while 
the non-relativistic treatment presents a universal behavior in the units of Figs. [7] In the limit of low densities the 
relativistic curves approach the universal non-relativistic curve. 

2) The relativistic treatment leads to values of the electron Fermi energy consistently smaller than the ones of 
the non-relativistic treatment. 

3) While in the non-relativistic treatment the electron Fermi energy can reach, by compression, infinite values 
as Rws ~^ 0, in the relativistic treatment it reaches a perfectly finite value given by Eq. (|37[) attained when Rws 
coincides with the nuclear radius R c . 

The universality of the electron Fermi energy with respect to the number of protons N p has been obtained by 
expressing the Coulomb potential energy eV in terms of the function (f> given by Eq. (|4]), and by introducing the scale 
factor b given by (JB]). Accordingly, the radius of the Wigner-Seitz cell has been expressed in terms of the nucleus 

radius ([T8|) which is proportional to N p ' . 

It is clear then, from above considerations, the relativistic treatment of the Thomas-Fermi equation introduces 
significant differences from the current approximations in the literature: a) the uniform electron distribution [22j . 
b) the approximate perturbative solutions departing from the uniform distribution [l4| and c) the non-relativistic 
treatment [16]. We have recently applied these results of the relativistic Feynman, Metropolis, Teller treatment of a 
compressed atom to the study of white dwarfs and their consequences on the determination of their masses, radii and 
critical mass (l5j . 



V. APPLICATION TO NUCLEAR MATTER CORES OF STELLAR DIMENSIONS 

We turn now to nuclear matter cores of stellar dimensions of A ~ (mpianck/^ra) 3 ~ 10 57 or M core <~ Mq. 
Following the treatment presented in Popov et al.,[l3j], we use the existence of scaling laws and proceed to the 
ultra- relativistic limit of Eqs. (|19p. (|2"6"|) . (|27p. (|28[), For positive values of the electron Fermi energy Ef , we introduce 

the new function (f> = 4 1 / 3 (97r)~ 1 / 3 xA/a; and the new variable x = kx where k = (12/t:) 1 ^ 6 ^/aA" 1 , as well as the 
variable £ = x — x c in order to describe better the region around the core radius. 
Eq. ([27]) becomes 



(44) 
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FIG. 8: The electron Coulomb potential energies in units of the pion rest mass in a nuclear matter core of stellar dimensions 
with A ~ 10 57 or Mcore ~ Mq and R c ~ 10 6 cm, are plotted as a function of the dimensionless variable £, for different values 
of the electron Fermi energy also in units of the pion rest mass. The solid line corresponds to the case of null electron Fermi 
energy. By increasing the value of the electron Fermi energy the electron Coulomb potential energy depth is reduced. 



where </>(£) = (/>(£ + x c ) and the curvature term 2</>' (£)/(£ + x c ) has been neglected. 
The Coulomb potential energy is given by 



Et 



corresponding to the electric field 



and the electron number-density 



1/6 



eh 



(45) 



(46) 



n e (r) 



1 



From Eqs. (fT9|) and (|47| we than have that, for £ = — x c , cj)(~x c ) = 1. 



(47) 



In the core center we must have n e = 

In order to consider a compressed nuclear density core of stellar dimensions, we then introduce a Wigner-Seitz cell 
determining the outer boundary of the electron distribution which, in the new radial coordinate £ is characterized by 
£ ws . In view of the global charge neutrality of the system the electric field goes to zero at £ = £ ws . This implies, 
from Eq. (gBJ), <j)'(£, WS ) = 0. 

We now turn to the determination of the Fermi energy of the electrons in this compressed core. The function <p 
and its first derivative <f>' must be continuous at the surface £ = of the nuclear density core. 

This boundary- value problem can be solved analytically and indeed Eq. (j44]) has the first integral, 



with boundary conditions at £ = 0: 



2[0'(O] 2 = 

0(0) 
0'(O) 



<t>\o - mo + 3, 

4>H0 - ^ ws ), 



£<o, 
e>o, 



(48) 



WS\ 



1 $ 4 (0) - ^^ ws ) 



(49) 
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FIG. 9: Solutions of the ultra-relativistic Thomas-Fermi equation (|44[) for different values of the Wigner-Seitz cell radius Rws 
and correspondingly of the electron Fermi energy in units of the pion rest mass as in Fig. near the core surface. The solid 
line corresponds to the case of null electron Fermi energy. 
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FIG. 10: The electric field in units of the critical field for vacuum polarization E c = m^c 3 /(eh) is plotted as a function of the 
coordinate £, for different values of the electron Fermi energy in units of the pion mass. The solid line corresponds to the case 
of null electron Fermi energy. To an increase of the value of the electron Fermi energy it is found a reduction of the peak of 
the electric field. 

Having fullfilled the continuity condition we integrate Eq. (|48|) obtaining for £ < 



m = 1-3 

where the integration constant a has the value 



1 + 2~ 1/2 sinh(a- V3£) 



sinh(a) = v2 



ii+^ 4 (r A ) 

1 - 04(£WS) 



(50) 



(51) 
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FIG. 11: The Fermi energy of electrons in units of the pion rest mass is plotted for different Wigner-Seitz cell dimensions 
(i.e for different compressions) ( ws in the ultra-relativistic approximation . In the limit S; WS —¥ the electron Fermi energy 
approaches asymptotically the value (E^) max given by Eq. ([62}. 



In the interval < £ < £ , the field 0(£) is implicitly given by 

F (arceos^p, ±\ = - £™) f (52) 

where F(ip, k) is the elliptic function of the first kind, and F(0, k) = 0. For F(ip, k) = u, the inverse function 
ip = F _1 (u, k) = am(u, k) is the well known Jacobi amplitude. In terms of it, we can express the solution (|52p for 
£ > as, 



m = 0(r 6 ) cos 



am(^r s )(£-r s ),i 



(53) 



In the present case of > the ultra-relativistic approximation is indeed always valid up to £ = £, ws for high 
compression factors, i.e. for Rws — Rc- In the case E^ = 0, £ ws — > oo, there is a breakdown of the ultra-relativistic 
approximation when £ — > £ ws . 

Details are given in Figs. 13 [TO] 

We can now estimate two crucial quantities of the solutions: the Coulomb potential at the center of the configuration 
and the electric field at the surface of the core 

^(O^^Y 73 ^ 2 -^, (54) 



£ max ~ 2.4^^J E e \$(0)\, (55) 

where E c — m 2 e C" / (eh) is the critical electric field for vacuum polarization. These functions depend on the value 
<i>{£ ws ) via Eqs. (|4"5 |) -([52 ]) . At the boundary £ = £, ws , due to the global charge neutrality, both the electric field 
E(£, ws ) and the Coulomb potential eV(^ ws ) vanish. From Eq. (|4"5]l. we determine the value of <j>(£) at £ = £, ws 

0(r S ) = Af±) 1/3 ^, (56) 



17 



as a function of the electron Fermi energies E F . From the above Eq. (l56|) . one can see that there exists a solution, 
characterized by the value of electron Fermi energy 



(£f)r 



1 /9tt 

a It 



such that 0(£ w s ) = 1. From Eq. (JHSJ) and £ = 0, we also have 



1/3 



(57) 



0(0) 



arccos 4 



For 0(£ vvs ') = 1, from Eq. (gU) follows 0(0) = 1 hence Eq. 



^ S (0(O)) = F 



0. 



0(0) ( 

becomes 
1 

71 



i 

71 



It is well known that if the inverse Jacobi amplitude F[0, 1/V2] is zero, then 

£WS(fcW S) = ^ (0) = 1} = o. 



(58) 



(59) 



(60) 



Indeed from j){£ ws ) = 1 follows 0(0) = 1 and £ ws = 0. When £ ws = from Eq. (@9j follows 0'(O) = and, using 
Eq. ([55]) , E max = 0. In other words for the value of E F fulfilling Eq. (|56]l no electric field exists on the boundary of 
the core and from Eq. (|4"7} and Eqs. (|18[ IT9")) it follows that indeed this is the solution fulfilling both global N e = N p 
and local n e = n. p charge neutrality. In this special case, starting from Eq. (|28|) and A = N p + N n , we obtain 

^L(hc) 3 A- - (E F ) 3 

(J?F\3/2 _ 4 y ,lC ) W a )max 

K-^e /max ,. _ ,_ -.3/4' l ul J 



2 3/2 



) 3 



2/3 



+ mlc 4 



In the ultra-relativistic approximation (E F ) max /^-(fic) 3 ^ << 1 so Eq. (foTj) can be approximated to 

2/3 



F ) 

e Jmax 



2 l/3< 



-1 



'1 



2 7 3 



TOttC 2 , 



where 



i^' = V'1 



(62) 



(63) 



The corresponding limiting value to the N p /A ratio is obtained as follows 

i 1 2 



i4 



2 7 3 



-1 



'1 



2 7 3 



(64) 



Inserting Eqs. (|62l) . (jSBf in Eq. (|64l) one obtains the ultra-relativistic limit of Eq. (|37p . since the electron Fermi energy, 
in view of the scaling laws introduced in [l3j , is independent of the value of A and depends only on the density of the 
core. 

The iVp-independence in the limiting case of maximum electron Fermi energy attained when Rws — Rci in which the 
ultra-relativistic treatment approaches the uniform one, and the iVp-dependence for smaller compressions Rws > Rc 
can be understood as follows. Let see the solution to the ultra-relativistic equation (|4*4l for small £ > 0. Analogously to 
the Feynman-Metropolis- Teller approach to the non-relativistic Thomas- Fermi equation, we solve the ultra-relativistic 
equation (|44|) for small £. Expanding 0(£) about £ = in a semi convergent power series, 



0(0) 



= i + £; a„e 



/2 



(65) 
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and substituting it into the ultra- relativistic equation (|44p . we have 



exp 



fc=3 



31n(l + ^a„r /2 ) 



71=2 



(66) 



This leads to a recursive determination of the coefficients: 

a 3 = 0, a 4 = 2 (O)/2, a 5 = 0, a 6 = 2 (O)a 2 /2, a 7 = 0, a 8 = </> 2 (0)(l - al)/8, • • •, (67) 
with a2 = ft (0) / (f)(0) determined by the initial slop, namely, the boundary condition ft(0) and (f)(0) in Eq. (|49l) : 



6'(0) = -1 



(68) 



Thus the series solution (1551) is uniquely determined by the boundary value (f)(^ ws ) at the Wigner-Seitz cell radius. 
Now we consider the solution up to the leading orders 



m = m + ftm + l^me + l^m^e + l^ma - *i)e + 



Using Eq. (fBT)]) . the electron Fermi energy (|5B1) becomes 

E e = ( E e)rr, 



i + « 2 r s + \pma ws ? + ^ 2 (o)« 2 (r s ) 3 + l$ 2 (o)(i - aM ws r + 



(69) 



0(0), (70) 



where (Ef) max = (97r/4) 1 / 3 A 1 is the maximum Fermi energy which is attained when the Wigner-Seitz cell radius 
equals the nucleus radius R c (see Eq,l57). For $(£, ws ) < 1, we approximately have (f>(0) = 3/4, ft(0) = -(3/4) 2 /\/2 
and the initial slope a 2 = ft (0) / 4>(0) = -(3/4)/V2. Therefore Eq. d70j) becomes 



E 



±V2 



\vs 



y 2 ) 



(71) 



By the definition of the coordinate £, we know all terms except the first term in the square bracket depend on the 
values of N p . In the limit of maximum compression when the electron Fermi energy acquires its maximum value, 
namely when £ ws — 0, the electron Fermi energy (|71l) is the same as the one obtained from the uniform approximation 
which is independent of N p . For smaller compressions, namely for £ ws > the electron Fermi energy deviates from 
the one given by the uniform approximation becoming iV p -dependent. 

In Fig. [11] we plot the Fermi energy of electrons, in units of the pion rest mass, as a function of the dimensionless 
parameter ^ w/,s and, as £ ws — > 0, the limiting value given by Eq. (|62"T) is clearly displayed. 

In ref. [§), in order to study the electrodynamical properties of strange stars, the ultra-relativistic Thomas- Fermi 
equation was numerically solved in the case of bare strange stars as well as in the case of strange stars with a crust 
(see e.g. curves (a) and (b) in Fig. 6 of ref. [§]). In Fig. 6 of @ was plotted what they called the Coulomb potential 
energy, which we will denote as VAicock- The potential VAicock was plotted for different values of the electron Fermi 
momentum at the edge of the crust. Actually, such potential VAicock is not the Coulomb potential eV but it coincides 
with our function eV — eV + . Namely, the potential VAicock corresponds to the Coulomb potential shifted by the 
the Fermi energy of the electrons. We then have from Eq. (|43|) 



eV(0 



1/3 



A 



m^c 2 (f)(0 = VAicock. 



(72) 



This explains why in [9( , for different values of the Fermi momentum at the crust the depth of the potential VAicock 
remains unchanged. Instead, the correct behaviour of the Coulomb potential is quite different and, indeed, its depth 
decreases with increasing of compression as can be seen in Fig. |8] 



VI. COMPRESSIONAL ENERGY OF NUCLEAR MATTER CORES OF STELLAR DIMENSIONS 

We turn now to the compressional energy of these family of compressed nuclear matter cores of stellar dimensions 
each characterized by a different Fermi energy of the electrons. The kinematic energy-spectra of complete degenerate 
electrons, protons and neutrons are 



l ( P ) - yW+^c 4 , V<Pf, i = e,p,n. (73) 
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So the compressional energy of the system is given by 



£ -\~ £ G ~f~ ^em 5 

d 3 rd 3 p 



Si 



E 2 



-d A r. 



(27Tft) 3 

Using the analytic solution (|53p we calculate the energy difference between two systems, I and 77, 

A£ = £(EF(II))-£(eF(I)), 



(74) 
(75) 

(76) 



with El (II) > El (I) > 0, at fixed A and R c . 

We first consider the infinitesimal variation of the total energy 6£tot with respect to the infinitesimal variation of 
the electron Fermi energy 8El 



S£ 



For the first term of this relation we have 
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(77) 



(78) 



where the general definition of chemical potential dei/drii = d£i/dNi is used (i = e, p, n) neglecting the mass defect 
m n — nip — m e . Further using the condition of the beta-equilibrium (f2"5f we have 



(79) 
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For the second term of the Eq. ([77| we have 
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(80) 



since in the process of increasing the electron Fermi energy namely, by decreasing the radius of the Wigner-Seitz cell, 
the system by definition maintains the same number of baryons A and the same core radius R c . 
Now S£ reads 
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(81) 



so only the electromagnetic energy and the electron energy give non-null contributions. 
From this equation it follows that 



A£ = A£ c 



A£ c 



(82) 



where A£ cm = £ cm (El(II)) - £ cm (El(I)) and A£ c = £ e (El(II)) - £ c (El(I)). 
In the particular case in which El (II) = (El) max and El (I) = we obtain 



A£ ~ 0.75- 



3 5/3 /n l/3 l , % sl/Q 



p 



(83) 



which is positive. 

The compressional energy of a nuclear matter core of stellar dimensions increases with its electron Fermi energy as 
expected. 
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VII. CONCLUSIONS 

We have generalized to the relativistic regime the classic work of Feynman, Metropolis and Teller, solving a com- 
pressed atom by the Thomas-Fermi equation in a Wigner-Seitz cell. 

In the relativistic generalization the equation to be integrated is the relativistic Thomas- Fermi equation, also called 
the Vallarta- Rosen equation [2!| . The integration of this equation does not admit any regular solution for a point- like 
nucleus and both the nuclear radius and the nuclear composition have necessarily to be taken into account 0, [f| . 
This introduces a fundamental difference from the non-relativistic Thomas-Fermi model where a point-like nucleus 
was adopted. 

As in previous works [H, [f| fl2j the protons in the nuclei have been assumed to be at constant density, the electron 
distribution has been derived by the Thomas-Fermi relativistic equation and the neutron component has been derived 
by the beta equilibrium between neutrons, protons and electrons. 

We have also examined for completeness the relativistic generalization of the Thomas-Fermi-Dirac equation by 
taking into due account the exchange terms adopting the general approach of Migdal, Popov and Voskresenskii 
Q, and shown that these effects, generally small, can be neglected in the relativistic treatment. 

There are marked differences between the relativistic and the non-relativistic treatments. 

The first and the most general one is that the existence of a finite size nucleus introduces necessarily a limit to the 
compressibility: the dimension of the Wigner-Seitz cell can never be smaller then the nuclear size. Consequently the 
electron Fermi energy which in the non-relativistic approach can reach arbitrarily large values, reaches in the present 
case a perfectly finite value whose expression has been given in analytic form. There are in the literature many papers 
adopting a relativistic treatment for the electrons together with a point-like approximation for the nucleus, which is 
clearly inconsistent (see e.g. [H, H3]). 

The second is the clear difference of the electron distribution as a function of the radius and of the nuclear composi- 
tion as contrasted to the uniform approximation often adopted in the literature (see e.g.(22j) which we have explicitly 
shown in the Fig. [5] of Sec. IIVI Inferences based on the uniform approximation are not appropriate both in the 
relativistic and in the non-relativistic regime. 

The third, one of the most relevant, is that the relativistic Feynman-Metropolis- Teller treatment allows to treat 
globally and in generality the electrodynamical interaction within the atom and the relativistic corrections leading to 
a softening of the dependence of the electron Fermi energy on the compression factor, as well as a gradual decrease 
of the exchange terms in proceeding from the non-relativistic to the fully relativistic regimes. It is then possible to 
derive, as shown in Table[I]of Sec. lIVI a consistent equation of state for compressed matter which overcomes some of the 
difficulties of existing treatments describing the electrodynamical effect by a sequence of approximations which have 
lead to the occurrence of unphysical regimes e.g. the existence of negative pressure as in the Salpeter approach. As a 
direct application of this treatment we have reconsidered the study of white dwarfs within the relativistic Feynman, 
Metropolis, Teller approach and evaluate their effects on the value of the radii, of the masses of the equilibrium 
configurations as well as on the numerical value of the critical mass (l5| . We have there compared and contrasted the 
results obtained by Chandrasekhar with a uniform approximation with the ones obtained by the equation of state of 
Salpeter and the ones following from the treatment presented in this article. 

We have then extrapolated these results to the case of nuclear matter cores of stellar dimensions for A « 
(wpianck/?7in) 3 ~ 10 57 or M core ~ M & . The aim here is to explore the possibility of obtaining for these systems 
a self-consistent solution presenting global and not local charge neutrality. The results generalize the considerations 
presented in the previous article corresponding to a nuclear matter core of stellar dimensions with null Fermi energy 
of the electrons [13j . The ultra-relativistic approximation allows to obtain analytic expressions for the fields. The 
exchange terms can in this approximation be safely neglected. An entire family of configurations exist with values of 
the Fermi energy of the electrons ranging from zero to a maximum value (E^) max which is reached when the Wigner- 
Seitz cell coincides with the core radius. The configuration with = (E^) max corresponds to the configuration 
with N p — N e and n p — n e . For this limiting value of the Fermi energy the system fulfills both the global and the 
local charge neutrality and correspondingly no electrodynamical structure is present in the core. All the other config- 
urations presents overcritical electric fields close to their surface. The configuration with Ef = has the maximum 
value of the electric field at the core surface, well above the critical value E c (see Fig. [5J Fig. [S] and Fig. [TO] of Section 
|V)) . All these cores with overcritical electric fields are stable against the vacuum polarization process due to the Pauli 
blocking by the degenerate electrons [30(. We have also compared and contrasted our treatment of the relativistic 
Thomas- Fermi solutions to the corresponding one addressed in the framework of strange stars [9( pointing out in these 
treatments some inconsistency in the definition of the Coulomb potential. 

We have finally compared the compressional energy of configurations with selected values of the electron Fermi 
energy. In both systems of the compressed atoms and of the nuclear matter cores of stellar dimensions a maximum 
value of the Fermi energy has been reached corresponding to the case of Wigner-Seitz cell radius Rws coincident with 
the core radius R c . 
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Both problems considered, the one of a compressed atom and the one of compressed nuclear density core of stellar 
dimensions, have been treated by the solution of the relativistic Thomas- Fermi equation and by enforcing the condition 
of beta equilibrium. They are theoretically well defined and, in our opinion, a necessary step in order to approach the 
more complex problem of a neutron star core and its interface with the neutron star crust. 

Neutron stars are composed of two sharply different components: the liquid core at nuclear and/or supra- nuclear 
density consisting of neutrons, protons and electrons and a crust of degenerate electrons in a lattice of nuclei [H, H3] 
and possibly of free neutrons due to neutron drip when this process occurs (see e.g. [l|). Consequently, the boundary 
conditions for the electrons at the surface of the neutron star core will have generally a positive value of the electron 
Fermi energy in order to take into account the compressional effects of the neutron star crust on the core [44J . The 
case of zero electron Fermi energy corresponds to the limiting case of absence of the crust. 

In a set of interesting papers [3lM36| Glendenning has relaxed the local charge neutrality condition for the description 
of the mixed phases in hybrid stars. In such configurations the global charge neutrality condition, as opposed to the 
local one, is applied to the limited regions where mixed phases occur while in the pure phases the local charge 
neutrality condition still holds. We here generalize Glendenning's considerations by looking to a violation of the local 
charge neutrality condition on the entire configuration, still keeping its overall charge neutrality. This effect does not 
occur in processes occurring just locally, but it requires the global description of the equilibrium configuration. As 
exemplified in [HI, [43|, where gravitational effects are taken into account both in the T = and T 7^ cases, the 
local properties need the previous knowledge of the entire equilibrium configuration and therefore the description is 
necessarily global. 

In conclusion the analysis of compressed atoms following the relativistic Feynman, Metropolis, Teller treatment 
presented in the first part of this article has important consequences in the determination of the mass-radius relation 



of white dwarfs 15[ leading to the possibility of a direct confrontation of these results with observations, in view of the 
current great interest for the cosmological implications of the type la supernovae [38144 1| . The results presented in the 
second part of this article on nuclear matter cores of stellar dimensions evidence the possibility of having the existence 
of critical electromagnetic fields in the interface of the core and the neutron star crust. The results here obtained in 
a simplified but rigorous approach of the application of the relativistic Feynman, Metropolis, Teller treatment to the 
constant density cores in beta equilibrium conform to the prediction of Baym, Bethe and Pcthick 1]. This treatment 
has been further extended to the case in which a self- gravitating system of degenerate neutrons, protons and electrons 
is considered within the framework of relativistic quantum statistics and Einstein-Maxwell equations [44TJ and to the 
case in which also strong interactions are present [45[. 
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